Comment on the paper 'Inelastic granular gas: Computer simulations and kinetic 
theory of the cooling state' by J.M. Salazar and L. Brenig 



On 
On 

> 

o 
oo 



43 

o 



CO 

a 

C 

o 
o 



> 



On 
On 

c3 



C 

O 

o 



X 
S3 



Yves Elskensu 
Equipe turbulence plasma, case 321, 
Laboratoire de physique des interactions ioniques et moleculaires, UMR 6633 CNRS-Universite de Provence, 
Centre de Saint- Jerome, av. escadrille Normandie-Niemen, F-13397 Marseille Cedex 20 

(preprint TP99.14) 



A hierarchy of equations was introduced recently to de- 
scribe the kinetics of homogeneous cooling of the gas of in- 
elastic hard spheres. It is argued that this hierarchy does not 
describe this system accurately, and a simple test for this is 
proposed. A stochastic interpretation of the hierarchy is pre- 
sented. 



PACS numbers: 
83.70.Fn Granular solids 
05.20.Dd Kinetic theory 



In their recent papers |lj|2j], Brenig and Salazar inves- 
tigate the kinetic theory of inelastic hard spheres. They 
obtain a hierarchy of ordinary differential equations for 
the moments Xk = (W k ) = (|fi2 ■ vi2| fe ) of the modulus 
W of the radial component of the relative velocity of two 
particle in the homogeneous cooling regime. Their equa- 
tions are remarkably simple and may shed new light on 
this fundamental process of granular media. 

It seems however that their derivation is incorrect. 
This Comment aims at (1) introducing a probabilistic 
interpretation of their equations, (2) arguing on the core 
of the illegitimate derivation and (3) on a first-principle 
problem raised by it, (4) commenting on processes for 
which their equations may hold and (5) proposing a sim- 
ple test for these equations. 

1. The Brenig-Salazar hierarchy reads for all k > 1 



ik = -(1 - a k )xk-\ 



(1) 



with the dot denoting derivation with respect to rescaled 
H time t — t'/r, where t' is the physical time, and r is the 
characteristic intercollision distance. Given the spheres 
diameter a and number density n, r = (27ra 2 n) _1 in 
three space dimensions and t = (iran)^ 1 in two space 
dimensions. The dimensionless parameter < a < 1 is 
the (normal velocity) restitution factor at the collision. 

Denoting by f(w,t) the probability distribution func- 
tion of W, hierarchy (|l|) is meaningful provided all mo- 
ments Xkit) — (W k ) — J °° w k f(w,t)dw are finite. Then 
f(w,t) itself is uniquely determined tf by the moments 
Xk{t). Elementary algebra rewrites (JjJ) as 



dtf(w,t) — —wf(w,t) + a 2 wf(a w,t) 



(2) 



The latter equation admits a natural interpretation as a 
random walk over [0, +oo[ for the normal relative velocity 
W of a pair of particles : for At — > 0, if W(t) = w, it may 
jump to W(t + At) — aw with probability wAt during 
the time interval At, and it may remain W(t + At) = w 
with probability 1 — wAt. Thus successive velocity jumps 
(collisions) occur independently, with a waiting time be- 
tween them distributed according to an exponential law 
(with parameter W(t)). 

2. Following notations of the appendix of Ref. g], this 
stochastic process describes indeed the result of inelastic 
collisions between particles 1 and 2. 

However both particles are surrounded by other par- 
ticles, and a collision of (say) particle 1 with a third 
particle 3 will change V12 to some v' 12 : in general 
I?i2 ' 1 7^ |?i2 • v^ 2 1 . Salazar and Brenig claim that 
the balance of these processes vanishes. This is doubful 
mathematically and physically. 

The mathematical difficulty originates from the Heav- 
iside functions 9(±gi3) in (A7) and from the fact that 
Iffial ^ I512I, l^ial ^ \9iz\- The two terms in the integrand 
cannot be reduced to each other, thus only a deeper sym- 
metry can yield = 0. The alleged symmetry, that 
expression (A10) be odd in fi3, overlooks the collisional 
requirement that g' 13 < || : (A10) contributes to (A7) 
only for a half-space integration, with constant sign. The 
other half-space integration also contributes with a con- 
stant sign, with a different integrand. Nothing warrants 
a cancellation. 

Physically, for a = 1, in a homogeneous stationary 
state, detailed collisional balance ensures that processes 
13 — * 1'3' have rates equal to processes 1'3' — > 13. But 
such a detailed balance does not apply to inelastic colli- 
sions - even in a homogeneous cooling regime. || 

3. From a more general viewpoint, the probability dis- 
tribution function of v' 12 depends not only on velocities 
vi , V2 and V3 but also on the distribution of the radial 
vector at impact fi3. Assuming that the distribution of 
the latter is continuous, the projection W — |f 12 • v' 12 | 
after collision 13 has a continuous probability distribu- 
tion ; one must even allow for W > W since particle 3 
may hit particle 1 on any side. Third-party collisions thus 
contribute to dtf an integral kernel which is continuous 
(rather than the Dirac distribution S(w- ln — a _1 ui ou t) of 
(^)). Now, the right hand side of (||) is a sum of terms de- 
scribing the various collisions, each with positive kernel. 



1 



Once a process has been allowed, no competing process 
can compensate it statistically : averaging with other 
processes acts in a convex manner on the infinitesimal 
generators of the random process. 

Equation (^) accounts neither for processes increasing 
W nor for processes smoothing a discrete initial distribu- 
tion of W. Therefore is appears unable to describe accu- 
rately the kinetics of the gas of inelastic hard spheres. 

Adding to the right hand side of @ a diffusion term ||] 
seems a promising way to include the effect of third party 
collisions, though an integral kernel is more general. Note 
that, for TV 3> 1, collisions 13 and 23 are O(N) more fre- 
quent than direct collisions 12, so that their contribution 
would be diffusion-like by a central limit effect and could 
overwhelm the 12 term taken into account in (||). 

4. Does hierarchy |l]) model the inelastic Lorentz gas ? 
In this case, an ensemble of independent light particles 
move ballistically between collisions with fixed, infinitely 
heavy scattering obstacles ; the normal component of the 
light particle's momentum is multiplied by a at each col- 
lision. Again, the outgoing normal velocity of the light 
particle after one collision is generally not its incoming 
normal velocity for the next collision, as successive nor- 
mal vectors are not parallel, hence equations (|l|) and (||) 
do not apply. Moreover, with fixed scatterers of given 
size, the light particle may undergo a sequence of cor- 
related collisions, while equations (Q) and (||) describe 
independent collisions. 

There may however be a limit in which the differential 
equation (||) and hierarchy (|l|) could describe the evolu- 
tion of a light particle velocity. This limit is a memory- 
less stochastic model, in which scatterers are Poisson dis- 
tributed in space, the modulus (instead of normal com- 
ponent) of the light particle outgoing velocity is a times 
its incoming velocity modulus, and the outgoing direc- 
tion is chosen so that recollisions would have vanishing 
probability in a Grad limit (scatterer density n — > oo, 
scatterer size a — > 0, and t constant). This model is not 
a limit of the previous Lorentz gas. 

5. Finally, independently from the derivation of the hi- 
erarchy from the iV-body Liouville equation |^] , one may 
also test the relevance of hierarchy |l]) and equation (|^) 
to homogeneous cooling (or to other physical processes) 
directly through an explicit prediction. For kinetics de- 
scribed by (H), all expectations 



invitation to universite libre de Bruxelles which initiated 
these discussions. 



Cn(t) = 
Sn(t) = 



cos(2imh\{Kw)/\n.a)f{w,t)dw (3) 
sin(2Tmln(Kw)/ hi a) f(w,t)dw (4) 



are constant in time, where an arbitrary constant k is 
introduced for dimensional reasons. Indeed, the random 
walk associated with (||) changes \n(nW(t))/ In a by unit 
steps +1, leaving the trigonometric functions invariant. 

It is a pleasure to thank L. Brenig, R. Lambiotte and 
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